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CHANNEL NORM-BASED ORDERING AND WHITENED DECODING 

For 

MIMO COMMUNICATION SYSTEMS 

5 Related Applications 

This application is related to provisional application serial number 60/404,860, 
docket number TI-35082PS, filed 08/21/02, by David J. Love, Srinath Hosur, Anuj Batra, 
Tarik Muharemovic and Eko N. Onggosanusi. 

10 Background of the Invention 

1. Field of the Invention 

This invention relates generally to multiple-input multiple-output (MIMO) 
communication systems. More particularly, this invention relates to channel norm-based 
15 ordering and whitened decoding for MIMO communication systems. 

2. Description of the Prior Art 

Multiple-input multiple-output (MIMO) communication systems provide gain in 
capacity and quality compared to single-input single-output (SISO) communication 

20 systems. While a SISO system employs one transmitter and one receiver to receive the 
transmitted signal, a MIMO system in general uses M t >= 1 transmitters and M r >=1 
receivers. Thus the SISO system is a special case of a MIMO system, with M t = M r =1 . 
Examples of MIMO systems include but are not limited to: 1. A communication system 
employing multiple-antennas at the transmitter and/or receiver; 2. A communication 

25 system employing orthogonal frequency division multiplexing (OFDM) or code division 
multiplexing (CDMA); 3. A time/frequency division multiple access system; 4. Any 
multiuser communication system; 5. Any combination of 1-4 above. 

Typically, the MIMO systems consist of a MIMO transmitter that sends 
30 "multidimensional" symbol information. This multidimensional symbol could, but is not 
limited to, be represented by a vector (note that a matrix symbol can always be 
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represented as a vector symbol by stacking the columns/rows of the matrix into a vector). 
The multidimensional symbol might represent one or more coded or uncoded data 
symbols corresponding to SISO transmitters. The transmitted signal propagates through 
the channel and is received and processed by a MIMO receiver. Note that the receiver 
could obtain multiple received signals corresponding to each transmitted symbol. The 
performance of the communication system hinges on the ability of the receiver to process 
and find reliable estimates of the transmitted symbol based on the received signals. 

Definitions 

As used herein, bolded capital symbols, such as H, represent matrices. 
As used herein, bolded lower-case symbols, such as s, represent vectors. 
As used herein, T denotes matrix transposition. 
As used herein, * denotes the matrix conjugate transpose operation. 
As used herein, _/ denotes the matrix inverse operation. 
As used herein, if W is a matrix, W m denotes the mth column of W. 
As used herein, if W is a matrix, (W r ) m denotes the mth row of W. 
As used herein, if v is a vector, ||v|| 2 , denotes the 2-norm of v. 
As used herein, if QO represents the symbol slicing function, it will be assumed 
to slice both single symbols and multi-dimensional symbol vectors. 
As used herein, Im represents the MbyM identity matrix. 
As used herein, 0 Mx n represents the M by N matrix of zeros. 
As used herein, if A and B are sets, then A\B is the set of all elements in A that are 

not in B. 

For MIMO systems such as, but not limited to, the ones discussed herein above, 
the received signal can be written, after front end receive processing such as filtering, 
downconversion, AGC, synchronization etc., in the form 
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where H„ is an M r by M t matrix of complex gains, s* is the M t -dimensional symbol 
vector transmitted at time k, and v is a M t -dimensional vector of additive noise. In 
narrowband wireless systems where the symbol period is much larger than the RMS 
delay spread as well as in OFDM systems where the inter-symbol interference is 
5 negligible due to the insertion of a cyclic prefix and/or guard interval, the channel from 
each transmit antenna to each receive antenna (per frequency bin in case of OFDM) is 
often modeled as a single-tap complex gain. In this case equation (1) simplifies to 



10 



y* = Hs* + v ^ 2 ) 



where H is now an M r by M t matrix of complex numbers and Hs* is the matrix product of 
H and s*. 



The receiver must estimate the symbol matrix S=[si ... s^] in order to facilitate 
15 reliable communication. Examples, but by no means the only examples, of 

multidimensional symbols could be space-time codes where 7>1 or spatial multiplexing 
systems with T=\ and independent SISO modulation on each transmit antenna. In case of 
no or negligible additive noise, v, and an invertible H, the estimation problem would 
reduce to that of inverting H. The presence of non-negligible noise, however, increases 
20 the difficulty in estimating S. Note that we have assumed that the receiver is has some 
estimate of H, that could be obtained by transmitting appropriate training sequences. The 
symbol matrix S is also assumed to be chosen from a finite set C of possible 
multidimensional symbols (this is typically the case as for e.g. where each element of S is 
chosen from a QAM symbol set). 

25 

The optimal solution in the sense of minimizing the probability of symbol error 
has been shown to be the maximum aposterior (MAP) decoder which in case of 
equiprobable symbol transmissions is equivalent to a maximum likelihood (ML) decoder. 

The ML decoder attempts to find S, the symbol matrix, by using the symbol matrix S 

30 that maximizes /?^S|yi,....,yrj where p (. |y i y *r ) is the conditional probability 
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density function (pdf) of s* given yi,..., y r . In real-time communications systems, 
however, this type of decoder is overly computationally complex. Decoders that search 
over a set V of possible multidimensional symbols S and decode to the multidimensional 

symbol S in Fthat minimizes some sort of metric are denoted as minimum distance 
(MD) decoders. The MAP and ML decoders are MD decoders with V=C 9 where C is the 
set of all possible multidimensional symbols S. 

Many algorithms that are computationally easier than ML decoding have been 
proposed in order to overcome the huge computational burden of ML decodings. 
Algorithms that perform some form of reduced complexity decoding will be referred to 
herein as sub-optimal decoders. An example of a suboptimal decoder is successive 
interference cancellation (SIC) method. A receiver using SIC decodes each symbol 
within the symbol vector one at a time. After each symbol is decoded, its approximate 
contribution to the received vector is subtracted in order to improve the estimate of the 
next symbol within the symbol vector. The order of symbol decoding and subtraction 
could be arbitrary or based on rules such as maximization of pre/post processing SNR 
etc. 

An example of an SIC receiver is the ordered iterative minimum mean squared 
error (IMMSE) receiver. With a single-tap channel, the receive signal is given by 
equation (2) above. Letting s k = |>i s 2 ... s Mr ] T , the ordered IMMSE operates using the 
following steps, letting y>,o = y*,Z>o = {l,2,..., M}, and H* (0) = H . 
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1. Setm = 0. 



2. Compute W m = (Hj*H 4 w + phf, - m )' l H^ . 
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6. Repeat steps 1-5 for m<M t . 

7. Set the decoded symbol vector to s* = sk,2 ... 7 mm if . 

Regarding the above algorithm, it is important to note that H^ +1) denotes the ith row of 

5 the matrix H* m+1) (time k and iteration m+i). Another example among sub-optimal 

decoders is the zero-forcing decoder which decodes to the symbol sk = £?(H -1 y*). This 
decoder is usually considered the worst performing and least complex of the sub- 
optimal decoders. 

10 Sub-optimal techniques unfortunately differ in diversity order from ML decoding 

(i.e. the asymptotic slope of the average probability of bit error curve). They essentially 
trade reduced complexity for reduced performance. 

In view of the foregoing, it is both advantageous and desirable to provide a lower 
15 complexity iterative decoder with channel norm based ordering that performs 

approximately the same as the iterative minimum mean squared error (MMSE) decoder. 
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Summary of the Invention 

The present invention is directed to a channel norm-based ordering and whitened 
decoding technique (lower complexity iterative decoder) for MIMO communication 
5 systems and that performs approximately the same as the iterative minimum mean 
squared error decoder. 

According to one embodiment, a method of decoding a signal vector comprises 
the steps of receiving a signal vector y* ; multiplying the received signal vector y* by a 

10 conjugate transpose of a channel matrix H* and generating a column vector z* 

therefrom; reordering entries associated with the column vector z* and generating an 
estimated channel matrix H therefrom; decomposing the estimated channel matrix H 
via Cholesky decomposition and generating a triangular matrix L therefrom; solving 
triangular matrix L backwards and estimating a signal vector s* therefrom, wherein s* 

15 is the true sorted symbol vector; and sorting signal vector s* and generating an estimate 
of the transmitted symbol vector s* therefrom. 
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Brief Description of the Drawings 



Other aspects and features of the present invention and many of the attendant 
advantages of the present invention will be readily appreciated as the same become better 
5 understood by reference to the following detailed description when considered in 
connection with the accompanying drawings wherein: 

Figure 1 is a system block diagram illustrating a decoder according to one 
embodiment of the present invention; 



10 



Figure 2 is a graph depicting simulation results illustrating a comparison of the 
decoder shown in Figure 1 with an ordered IMMSE for a 2x2 system using 16-QAM; and 



Figure 3 is a graph depicting simulation results illustrating a comparison of the 
15 decoder shown in Figure 1 with an ordered MMSE for a 2x2 system using 64-QAM. 



While the above-identified drawing figures set forth particular embodiments, 
other embodiments of the present invention are also contemplated, as noted in the 
discussion. In all cases, this disclosure presents illustrated embodiments of the present 
20 invention by way of representation and not limitation. Numerous other modifications and 
embodiments can be devised by those skilled in the art which fall within the scope and 
spirit of the principles of this invention. 
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Detailed Description of the Preferred Embodiments 



Figure 1 is a system block diagram illustrating a decoding algorithm 100 
according to one embodiment of the present invention. The decoding algorithm 100 was 

5 found by the present inventors to provide reduced complexity without much loss in 
performance when compared with known sub-optimal schemes based on iterative 
interference cancellation. In order to preserve clarity and brevity, and to better 
understand how decoding algorithm 100 functions, it will be assumed that each random 
entry in v is complex Gaussian with variance a 2 and that the expected value of the 

10 transmit power on each antenna is constrained to be 1 . A received symbol vector y* is 
multiplied by the conjugate transpose of the channel matrix yielding 



15 The entries in the column vector z k are then reordered based on the column norms of the 
channel. This is done by defining 



z* = H*Hs* + H*v. 



(3) 




(4) 



20 and 



= [v*,/„v*./ J v t . <w J andfi = [ Ht . m . g^J 



(5) 



where for m < n , ||H*,*J 2 < |Hm„[| 2 . At this point, a Cholesky decomposition is taken 



25 such that 



L*L = H*H + cr 2 Lv 



(6) 



Now multiplying z* by L* 1 yields 



30 
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= L* '(iTS + a 2 lM^k - L^VIa/, S* + L*~ l H*\ = Ls* + L M (fi # v - (T 2 lM t Sk) (7) 

where s* is the true sorted symbol vector. The triangular matrix L is then used to solve 
backwards for s* . After estimating Sk , the vector is then sorted into an estimate of the 
5 transmitted symbol vector s* . The decoding algorithm 100 returns this estimated vector 
along with the soft and hard decoded bits as shown in Figure 1. 

It is important to note that the noise+interference term during the backwards 
substitution step, L*~ ! (h*v - <r 2 lM t Sk) 9 is zero mean with a covariance matrix of <j 2 lM t 
10 under the assumption that the transmitters transmit each point in the constellation with 
equal probability. This is important if an outer code such as a convolutional code is used 
since it allows the Viterbi algorithm to be used as an optimal ML decoding scheme for 
the outer code. The IMMSE algorithm does not have this property. 

15 The computational complexity of the decoding algorithm 100 shown in Figure 1 

is also lower than ordered IMMSE. This algorithm 100 requires ordering only once 
during the decoding. In contrast, the ordered IMMSE algorithm requires M f orderings 
while decoding a symbol vector. The algorithm 100 shown in Figure 1 and discussed 
herein above, can be seen to also require a Cholesky decomposition and the inversion of a 

20 triangular matrix. In contrast, the ordered IMMSE algorithm requires M f - 1 matrix 
inversions. 

The algorithm 100 was simulated for a system with 2 inputs and 2 outputs 
experiencing MEMO Rayleigh fading using quadrature amplitude modulation (QAM) for 

25 various constellation sizes. In MIMO Rayleigh fading, the channel matrix H is a random 
matrix where all the entries of H are independent, identically distributed complex 
Gaussian random variables. Simulations shown in Figure 2 and Figure 3 show that the 
raw (or uncoded) average probability of bit error of the algorithm 100 compares 
favorably with ordered IMMSE decoding. In fact, for 16-QAM 200, algorithm 100 can 

30 be seen to outperform ordered IMMSE decoding. 
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In summary explanation of the above, a column norm-based ordered, whitened 
decoding algorithm 100 is provided for MIMO communications systems. The decoding 
algorithm 100 has been shown to perform approximately the same as ordered iterative 
minimum mean squared error decoding while incurring less computational complexity. 
5 The decoding algorithm 100 as well does not color the noise between symbol decisions 
when decoding a symbol vector. 

In view of the above, it can be seen the present invention presents a significant 
advancement in the art of multiple-input multiple-output (MIMO) communication 
10 systems. Further, this invention has been described in considerable detail in order to 
provide those skilled in the decoding art with the information needed to apply the novel 
principles and to construct and use such specialized components as are required. 

Further, it should be apparent that the present invention represents a significant 
15 departure from the prior art in construction and operation. However, while particular 
embodiments of the present invention have been described herein in detail, it is to be 
understood that various alterations, modifications and substitutions can be made therein 
without departing in any way from the spirit and scope of the present invention, as 
defined in the claims which follow. For example, the although the inventive 
20 embodiments described herein before were described in association with a 2x2 system, 
the invention is not so limited, and can just as easily be applied to any MtxMr system. 
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